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Abstract 

Using a Ginzburg-Landau model, we study the vortex behavior of a rectangular 
thin film superconductor subjected to an applied current fed into a portion of the 
sides and an applied magnetic field directed orthogonal to the film. Through a center 
manifold reduction we develop a rigorous bifurcation theory for the appearance of 
periodic solutions in certain parameter regimes near the normal state. The leading 
order dynamics yield in particular a motion law for kinematic vortices moving up 
and down the center line of the sample. We also present computations that reveal 
the co-existence and periodic evolution of kinematic and magnetic vortices. 
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1 Introduction 

We consider a thin superconducting sample occupying a rectangle. The sample is sub- 
jected to normal electric current that enters through a lead on one side and leaves through 
another lead at the opposite side. In addition the sample is subjected to a magnetic field 
oriented in the direction perpendicular to the rectangle's plane. Our main interest is 
in the regime in the parameter space where different physical quantities, such as the 
total current in the sample and the order parameter, are time-periodic. We pay special 
attention to the formation and motion of vortices in the sample. 

The problem above is a natural generalization of the simpler case of a finite super- 
conducting one-dimensional wire subjected to normal current that is fed into one of its 
endpoints. This problem received considerable attention since it is a canonical case of 
co-existence of normal current and superconducting current. Moreover, it is known that 
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in this setting there exists a regime of prescribed current I and temperature T, where 
the sample's behavior is time-periodic. In addition, in this regime the superconducting 
order parameter ip(x,t) vanishes at the wire's center at specific points in time that are 
separated by a fixed period. Such zeros of ip are called phase slip centers (PSC). These 
phenomena and others are studied, numerically and experimentally, by a number of au- 
thors including [13, lHHS], and many of the results are summarized by Ivlev and Kopnin 

A recent study in [3T] and [52] presents a comprehensive theory that explains the 
different patterns observed in this wire setting. The key idea is that the underlying system 
of equations enjoys a PT symmetry, that is, symmetry under complex conjugation and 
the transformation x — > — x. This symmetry enabled the authors to perform a rigorous 
bifurcation study of the problem, and to deduce that, under certain conditions, the order 
parameter bifurcates, as the temperature is lowered beyond a critical value, from the 
normal state ip = to a nontrivial state. Moreover, this bifurcation is shown to be 
of Hopf type, and this explains the periodic nature of the solution and the periodic 
appearance of isolated zeros of ijj. 

In the present study we look at a more realistic geometry of a finite strip. Moreover, 
we consider not just forced electric current, but also the effect of an external magnetic 
field. The problem is analyzed numerically in [BJ. The authors observe, just as in the 
one-dimensional setting, a periodic behavior of a number of physical quantities, including 
periodic appearance and motion of vortices. Therefore, our goal is to derive a theory 
that explains the observed patterns and vortex motion. 

The issue of vortices is of particular interest. They are defined as isolated zeros of 
the order parameter, and they are characterized by their topological degree in the (x,y) 
plane. The appearance of vortices in superconducting samples subjected to an applied 
magnetic field is of course well-known. Therefore, we expect them in our setting even 
in the absence of forced electric currents. What makes the present problem interesting 
is that also the opposite is true; namely, vortices form, for appropriate range of values 
of I and T, even if no magnetic field is applied. Therefore, one can classify the vortices 
here into magnetic vortices generated by the applied magnetic field in the absence of any 
applied current, and kinematic vortices generated, as will be shown below, by the forced 
electric field and the special symmetry of the problem in the absence of any magnetic 
field. 

When both applied magnetic and electric fields are present, it is not as clear how to 
make a distinction between the two kinds of vortices. As will be shown below, the sym- 
metry of the problem implies that some vortices are formed and move time-periodically 
on the center line of the rectangle for large enough I and for a range of applied mag- 
netic fields h. We term them kinematic vortices. In certain cases, as presented below, 
these kinematic vortices collide and move off the center line line. We term such vortices, 
"born" from kinematic vortices, 'kinematic' as well. 

Other vortex phenomena, not usually observed in more standard Ginzburg-Landau 
settings, include the time-periodic emergence of vortex pairs that are of opposite de- 
grees, which we term vortex/anti- vortex pairs: typically stable vortex configurations of 
Ginzburg-Landau vortices only involve vortices of the same degree but in this periodic 
setting that is not always the case. We also show that for some range of I and h vor- 
tices of the same degree move towards each other and collide before moving away from 
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each other. Again, this is a process that is atypical to more familiar Ginzburg-Landau 
settings. 

The analysis of the present problem follows to some extent the lines of the one- 
dimensional wire problem. Namely, we construct a proper framework that enables us to 
use the Center Manifold Theorem to study the bifurcation picture, and thereby establish 
the existence of a Hopf bifurcation. Moreover, in both cases a key factor is played by 
the spectrum of the underlying linear Schrodinger operator. However, there are a few 
important differences between the one-dimensional case and the problem considered here. 
First, the construction of the center manifold requires certain apriori estimates on the 
solutions to the underlying differential equations. These estimates are harder to obtain in 
the two-dimensional setting. A second important difference relates to the vortex motion. 
While in the one-dimensional case the PSCs just appear momentarily at a fixed periodic 
sequence of instances, the kinematic vortices in the two-dimensional problem are present 
for periodic finite time inte rvals. Moreover, they tend to move along or near the y-axis 
(which is the center line of the rectangle) as we shall show. 

It is interesting to note that implications of PT-symmetry seem to appear in a number 
of quite different physical problems. For instance, we mention applications to quantum 
mechanics [4] , [8] , to hydrodynamic instability [24] , [25] , and to optics [19] , [20] . We also 
mention several recent rigorous studies within Ginzburg-Landau theory that incorporate 
magnetic effects along with applied currents in a variety of asymptotic regimes, [TJ [21 El 
[2S1I2S]- One aspect of our investigation that distinguishes it from others, however, is that 
it is to our knowledge the first to capture a motion law for Ginzburg-Landau vortices 
that is not based on the assumption of large Ginzburg-Landau parameter. 

In the next section we formulate the problem and the underlying equations. The bi- 
furcation analysis is performed in section 3. In particular we establish there the existence 
of a center manifold for a certain regime in the (J, T) plane. In section 4 we consider 
the formation of vortices and their motion and discuss some computational work on the 
problem. 

2 Formulation of Problem 

We consider a superconducting material occupying a thin rectangular box with dimen- 
sions — L < x < L, —K < y < K and say < z < r\ where rj is assumed to be much 
smaller than the coherence length or penetration depth, allowing us to work within the 
thin film 2d approximation of Ginzburg-Landau. In this approximation, we take the 
complex- valued order parameter = ^(x, y,t) and the real- valued electric potential 
4> = (j)(x, y, t) to be defined onKx [0, oo) where we denote 1Z := [-L, L] x {—K, K] and 
we ignore any induced magnetic field. Within this rectangular geometry, we assume the 
presence of leads forcing in electric current of magnitude / through the sides x = ±L, 
along the subinterval ~S < y < 6 for some positive 6 < K. Additionally we assume that 
the thin film is subjected to an applied magnetic field of size h oriented perpendicular to 
the rectangular cross-section. See Figure [T] 

Introducing the applied magnetic potential A := (— y, 0), this applied field is then 
given by hV x Aq and the thin film limit of Ginzburg-Landau takes the nondimension- 
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Figure 1: A thin film superconductor subjected to applied current and magnetic field, 
alized form 

^ t + i^= (V-i/i^ ) 2 * + (r-|*| 2 )* for (x, y) e U, t > 0, (2.1) 
Ac/) = V • f|{Wtf* - - \^\ 2 hA ^j for (x,y) eft, t > 0, (2.2) 

subject to the boundary conditions 

9(±L,y,t) = 0tor \y\ < S, (2.3) 
V x {±L,y, t) + thy*{±L,y, t) = for S <\y\^K, (2.4) 
$ y (x, ±K, t) = for |z| < L, (2.5) 

M±L,y, t) = {-' !0 %[>\<*- K (2.6) 
<j) y {x,±K,t) = for |ac| < L. (2.7) 

along with the initial condition ty(x,y, 0) = ^mit{x,y). As a convenient normalization 
we take J K = 0. The system of equations above are collectively known as the Time- 
Dependent Ginzburg-Landau (TDGL) model. 

The parameter T is proportional to T — T c where T is temperature and T c is the 
critical temperature below which the normal (zero) state loses stability in the absence 
of any applied fields. Note that (|2.2j) is simply the requirement of conservation of total 
(normal + superconducting) current. Also, we remark that the boundary conditions 
(12.4p - (|2.5p are the standard superconductor/vacuum conditions on "J, while (12.31) reflects 
the presence of the normal leads. 
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There are numerous investigations of current-driven superconducting wires and thin 
films that utilize a model based on Ginzburg-Landau theory, and these have generally 
reported reasonable agreement between theory and experiment. Regarding (|2.1I) . we 
should mention that some studies on this problem such as [SJ |H1 IE] replace the left-hand 
side with the modification 

d _ 7 2 3|*| 2 \ 



where u and 7 are material parameters. However, we believe that the standard and 
simpler evolution equation (I2.1[) corresponding to the choices u = 1 and 7 = captures 
the main features of the problem, an opinion shared by the authors of [3] whose compu- 
tational comparisons suggest that the modification does not have a large effect in this 
setting. What is more, in what follows, we shall concentrate on bifurcation from the 
normal state W = 0, so the smallness of the amplitude should mute the effect of this 
modification even more. 



3 Analysis of the model 

For a given "J, and a given value h 0, let us decompose the solution 4> to (|2.2j) . (|2.6[) . 
(12771) as 

4> = 14P + 4> 

and where <jP is harmonic and satisfies the boundary conditions 

^(±L,y) = h 1 , fo / ^<*V (3-1) 
xy \ for S < \y\ < K, y ' 

<j)° y (x, ±K) = for \x\ < L, (3.2) 

while <j> satisfies (|2.2|l subject to homogeneous Neumann boundary conditions on all 
portions of the rectangular boundary. We note that <j) but depends on ^ so we will often 
write <j) as 4>[^] to emphasize this dependence. 
We will occasionally make use of the properties 

cf> (-x,y) = -4>°(x,y) and <jP(x, -y) = <P°(x, y), (3.3) 

which are easy to check. 

The normal state in this setting corresponds to = and (j> = I<f>°. We will pursue 
a bifurcation analysis about this normal state, and therefore a crucial role will be played 
by the linear eigenvalue problem 

C{u] := (V - ihA ) 2 u - UcjPu = -Am for |ac| < L, \y\ < K, (3.4) 

subject to the boundary conditions f)2.3j) — (|2-5[) . Though we do not indicate it in our 
notation, it is understood that C and therefore all of its eigenvalues and eigenfunctions 
depend on the parameters L,K,6,h and /. We summarize below the key properties of 
the corresponding eigenvalues and eigenfunctions that will be needed in the analysis to 
follow. 
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Lemma 3.1. The spectrum of C consists only of point spectrum, denoted by with 
corresponding eigenf unctions {uj}. If(Xj,Uj) is an eigenpair satisfying (13.4j) then 

ReAj>0, and |Im A 3 -| < \\4> \\ L „ I- (3.5) 

Thus, in particular we may order the eigenvalues Ai, A2, • • ■ according to the size of their 

real part, with < Re Ai Re A2 ^ The PT-symmetry of the operator is reflected in 

the fact that if (Xj,Uj) is an eigenpair then so is (A*,itj) where Uj(x,y) := Uj(—x,y). 

When working with Miand U2, we will choose the normalization 

t ? = 1 = [ uj. (3.6) 
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Also, as a matter of convention, when Ai is non-real, we associate u\ with the eigenvalue 
Ai having positive imaginary part, and U2 with A2. 

Proof. The fact that the spectrum consists solely of eigenvalues follows from standard 
compact operator theory. The conditions in (|3.5p follow from multiplication of the equa- 
tion C[u] = —Xu by u* and integration over the rectangle. This leads to the identity 

In M In M 

implying (|3 . 5|) . 

The final claim of the lemma follows by noting that C[vJ] = C[u]. □ 



A numerical analysis of the eigenvalue problem Q3.4p indicates that, fixing all other 
parameters (i.e. K,L,S and h) and varying only the current J, there exists a critical 
value I c depending on these other parameters, such that Ai is real for I ^ I c and Ai is 
non-real for I > I c . This complexification of Ai arises through a collision with another 
eigenvalue and the two eigenvalues emerge from this collision, that is for / > I c , as 
complex conjugate pairs. Examples of eigenvalue collisions are shown in Figures and 
(2b • In both cases the simulated geometry is L = l,K = 2/3,5 = 4/15. In Figure^ we 
took h = 0, while in Figure [2b we took h = 7.5. 

We note in passing that in light of (|3 . 3[) and (|3.7p . the evenness of \uj\ evidently 
implies reality of Xj so the passage to complex eigenvalues with increased values of I 
carries with it a certain symmetry breaking of the corresponding eigenf unction. This 
phenomenon has been established rigorously for the one-dimensional version of (|3.4[) - 
that is, for the problem depending only on x and with h = in |24[ I25j . The reality of 
the spectrum for / positive and sufficiently small in this two-dimensional setting should 
follow by the type of perturbation analysis to be found in (7], since for I — 0, the spectrum 
is clearly real, cf. p.5[) . However, since here we are interested in the regime where the 
first eigenvalue is complex, we do not pursue this point further. 

Regarding eigenvalue collisions, we also wish to note another new phenomenon for this 
two-dimensional problem with the incorporation of magnetic effects that is not observed 
for the one-dimensional problem. Our computations reveal that for certain large enough 
values of applied magnetic field h, as I increases through the regime < I < I c , the 
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Figure 2: The real part of the spectrum of C. The first 4 eigenvalues are drawn for the 
parameters L = 1, K = 2/3,5 — 1/6. In (a) h — 0, while in (b) h — 7.5 
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Figure 3: The real parts of the four leading eigenvalues for the parameter values L = 
1, K = 2/3, 5 = 4/15, h = 20, / = 25. Note how A3 passes through A2 and then collides 
with Ai at I = J c « 20. 

second and third (still real) eigenvalues pass through each other and it is ultimately what 
was originally labeled as the third eigenvalue that collides with Ai at I c . We say the third 
eigenvalue "passes through" the second rather than "collides" with it because as / varies 
through the point where A2 = A3, both eigenvalues remain simple and the corresponding 
eigenfunctions vary smoothly without incident. This for example is the scenario for the 
parameter values L = 1, K = 2/3, 6 = 4/15, h = 20, / = 25. See FigureEl We will return 
to this set of parameter values at the end of the article to discuss anomalous vortex 
behavior as well. 

At this point we fix any values of K, L,5,h and then pick / sufficiently large so that 
ImAx 7^ 0. We then consider our problem (|2.ip - (l2.7[) with T given by T = ReAi + e 
where e is small and positive. Note that given the temperature dependence of T, this 
corresponds to lowering the temperature just below the value where the normal states 
becomes linearly unstable. 

We now introduce the linear operator L\ [u] := C[u] + (Re Ai) u. Recalling Lemma [37X1 
in the scenario where Ai is not real, the corresponding first eigenfunctions, say u\ and 
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U2, satisfy 

u 2 = u\ with Ci{u\\ = —ilm\\U\, C\\u,^ = ilm Ai «a- (3-8) 
For later use, we also note that if we introduce the adjoint operator C\ satisfying 

/ d[u]v= [ £l[v]u, 
Jtz Jtz 

then we can readily identify C*[v] as simply the operator given by (V + HiAq) 2 v — il<fi°v + 
(Re Ai) v. Furthermore, the fact that <po is even in y (cf. t\'3.3\i ) reveals that the functions 
u*(x, y) := Uj(x, —y) for j = 1, 2 satisfy the equations 

= — ilmAi u*, £^[1*2] = ilmAi «2- 

Consequently, we see that 

— ilmAi / u\u2 = j C\[u\\u2 — I C\ [u2\u\ — zlm Ai / u\u2, 
Jtz Jtz Jtz Jtz 

with a similar relation holding between u\ and Ui. Hence, ^u\u2 = = J-^u^ui. 
Clearly, this orthogonality holds between any two eigenfunctions iij and Uk of the operator 
C\ corresponding to distinct eigenvalues, namely 

i*u k = for j i= k. (3.9) 

TZ 

Of course, in the special case of no applied magnetic field, i.e. h — 0, we have C = C* 
and u\ = u* and through (|3.3|) we see then that each eigenfunction is even in y. 

Let us now return to the J2.1[) - (|2.2[) with the choice T = Re Ai + e and re-express the 
system as a single nonlinear, non-local equation for the complex-valued order parameter 

f t = A[*] +e^+Af(^), (3.10) 

where 

7V r (*) := -|*| 2 *-i0[*]*. (3.11) 

When expressed in this form, the problem can be rigorously solved for e <^ 1 via 
a center manifold reduction. The analysis is similar to that carried out for the one- 
dimensional (thin wire) problem in [22], Prop. 6.8, so we will only mention the key steps 
and those parts of the calculation where there are changes. 

As regards the linear part of (|3 . 10[) . the key point is that the operator — C\ is sectorial, 
in light of (13.51) . cf. [TT]. Regarding estimates on the cubic, nonlocal nonlinearity Af, 
the analysis differs from that in ;22] in that for the one-dimensional problem it is easy 
to check that Af is a bounded map from H 1 to H 1 (cf. [35], Lemma 6.3), while in the 
present two-dimensional setting this is no longer true-it just barely misses. One way to 
overcome this obstacle is by viewing Af as a mapping from an interpolation space between 
H 1 and H 2 into L 2 . The details of pursuing this strategy can be found in section 3.4 
of |17j , where the authors execute a center manifold construction relevant to bifurcation 
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from the normal state without applied electric current. Alternatively, one can view M as 
a map from H 2 (1Z) into H 1 ^). We describe how to make the necessary estimates for 
this latter approach. 

Writing Af = Aft + Af 2 with A/l($) := - |*| 2 * and A/^*) := it is an easy 

application of Holder's inequality to make the estimate 



Wi(n\m ( n)<C\\W HH K) 



(3.12) 



To make a similar estimate on A/2, let us first write the PDE coming from (|2.2[) for </>(\E') 

as 

(3.13) 



A0 = V- j(f)-|*l A 



subject to homogeneous Neumann conditions and the normalization condition of mean 
zero, where we have introduced the notation j(*ff) :— ^{WV^* — \f r *V\I r }. Fixing any 
p G (1, 2) we may use Holder's inequality and Sobolev imbedding to make the estimate 



|Dj(*)| p <C/ [ \#\ p \D 2 V\ P + \DV 

TZ JTZ 



2p 



I* 



,2p 

W 1 . 2 p(TC) 



<c\\n% {n) 



Hence, \\Dj(p)\\ 



Lp(K] 



H 2 (TZ) 



It is also easy to estimate 



\\m 



\Lp(K) 



m 2 A 



<C||*II^ 



(TZ) 



so we conclude that 



V - (jC*) - |*| 2 A ) 



Lp(K) 



(TZ) 



Then appealing to the Calderon-Zygmund inequality (cf. [27], Chapter 2), equation 
(I3.13[) implies that 



W 2 'P(TZ) 



H 2 (TZ) 



From the Sobolev imbedding theorem and Holder's inequality it then follows easily that 
11-^2 (*) || m(TZ) ^ C II*II'h 2 (7?,)- Combining this last estimate with (|3.12p we arrive at the 
estimate on the nonlinearity: 



\mn\ HH TZ)^c\\nl H TZ)- 



(3.14) 



The upshot is that for each small e, one can construct a center manifold M. e as a graph 
v M> $(w, e) in H 2 (1Z) over the center subspace S := span {m, U2}, applying for example, 
the version of the Center Manifold Theorem to be found in [TU], Theorem 2.9. More 
precisely, there exist positive constants Sq and £q, such that for any e satisfying |e| ^ Eq 



one can define M £ :— {$(«,£) : v G <S, 
properties: 



\H 2 (TZ) 



< S }, and A4 S enjoys the following 
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(i) The center manifold is locally invariant for the flow p.lOj) in the sense that if |e| < £q 
and the initial data tpo lies on M E , then so does the solution ^ e to p.lOp so long as 
IIV' e ( - 7 ^)\\h 2 (k) sta ys sufficiently small. Hence, for such initial data, one can describe the 
resulting solution tp £ (t) — tp £ (-,t) through two maps af, a| : [0, oo) — > C via ip £ (t) — 

(t)u\ + a§(t)w2) e). Since u\ and 112 are fixed throughout, we will often write simply 

r{t) = $>{a\(i),a e 2 {t),e) (3.15) 

for the sake of brevity. 

(ii) The center manifold is PT-symmetric, i.e. ip G A4 £ => G A4 £ . Furthermore, if 
v — v^, then e) — &(v, e). Since 112 — u\, this implies that if a|(0) = af (0)*, and one 
solves (|3.10p subject to the PT-symmetric initial conditions ip e (0) = $(af (0), a|(0), e), 
then the resulting functions a\(t) and a|(0 describing the solution ip £ at any positive 
time t will remain complex conjugates. 

(iii) A4 contains all nearby bounded solutions of (|3 . 10[) . and in particular, it contains 
any nearby steady-state or time-periodic solutions. 

(iv) Through an appeal to p,14j) . the discrepancy between the center manifold and the 
center subspace can be expressed through the estimate 

||$(«, e) — HI#2(tc) ^i(lMI^(7Z) + HIMW)) (3-16) 

which holds for any pair (v, e) such that v G S with |M| ff2 ( K ) < ^0 an d | £ | < £ 0: where 
C\ is a positive constant independent of v and e. 

Armed with these properties of the center manifold, one can then fix any sufficiently 
small complex numbers a\(0) and af (0) as in (i) above, solve (|3.10[) and then use the 
projection n c onto the center subspace S to obtain a reduced system of O.D.E.'s gov- 
erning the evolution of a\ and a\. In light of (|3.9[) we can write an arbitrary function / 

as / = C\U\ + C2U2 + u 1 - where w 1 G (span {u*, u!}) • Hence the projection onto the 
center subspace of an arbitrary function / is given by 



n c (/) = [f^r- «i+ [f^- u 2 . 

Thus, substituting the reduction (13.151) into (I3.10P and projecting, we obtain the following 
system through the use of (|3. 



a\u\ + a e 2 U2 = (e — £ImAi)a:fiti + (e + iImAi)a|u2 
+n c ( AT (aim + a E 2 u 2 ) ) + \ n c ( jV($(af,a|,£)) ] - II C ( jV(afui + a|u 2 ) 



(3.17) 

where ' denotes a time derivative. 

Invoking (|3.16[) . one finds that the last expression above involving the difference of 
nonlinear terms is lower order so one is justified in initially ignoring it, solving the 
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resulting simplified system of O.D.E.'s and then arguing that the behavior of solutions 
persists for the full system (|3.17p . Again the details of this portion of the argument can 
be found in [35]. In an abuse of notation, we will persist in using the notation to 
denote the solution to the truncated system in which the last expression is dropped. 

Now using (|3.9|) . we integrate first against u\ and then against itj to arrive at the 
system: 

■s i -T \ \ e , / tc <A/"Kmi + a|u 2 ) 
a\ = (e — ilm \i)a\ H — - 1 (3.18) 

■ e ( _i_ • T ^ e , S_K < + «1"2) Q s 

a 2 = (e + ihsx\\)an H ? — . (3.19) 

In order to argue that (|3.18l) - (|3.19l) exhibits a Hopf bifurcation to a periodic state, 
we now restrict the flow to the PT symmetric portion of the center subspace, and hence, 
in light of item (ii) above, to the PT symmetric subset of the center manifold. This 
amounts to the restriction a 2 = (af )* and allows us to only work with (|3.18j) . At this 
juncture, we make a change of variables of the form 

a E := a\ + Cl (af) 3 + c 2 \a\ | 2 (of)* + c 3 (K)*) 3 , 

in order to convert the problem to its normal form, cf. |10j . chapter 3. Writing the cubic 
nonlinear term in (I3.18P as 



ni(a? ) 3 + n 2 \a\\ 2 (of)* + n 3 (K)*) 3 + n 4 |af | 2 a 



a tedious but direct calculation yields that with the choices c\ = 2Ii '^ 1 Ai i, c 2 = — 2imfxT *' 
and C3 = — Tf^V" *i the new variable a e satisfies the simpler differential equation 



a E = (e - HmAiK + n 4 |a e | 2 a e , (3.20) 
with the coefficient n 4 given by 

~ Iv. ( U 1 U * + 2 M - « J* K ((011 + <I>22)U\U\ + (f>l2U*U 2 ) 

TI4 = p . (3.21) 

Here we have introduced the notation <pij to denote the solution to the equation 
Aipij = V • (-^[uiVu* — u*Vui] — u*UjhAo) = for i,j = 1,2 

subject to homogeneous Neumann boundary conditions on dlZ and zero mean on 1Z, so 
that <p{a e u\ + (a e )*v,2) in the nonlocal contribution A/2 to the nonlinearity Af takes the 
form |^| 2 (0 n + 22 ) + (a e ) 2 21 + ((a £ )*) 2 12 . 

Provided Re n 4 < 0, it is then easy to check from (I3.20[) that the system undergoes a 
supercritical Hopf bifurcation to a periodic solution given by 

a e( f ) = J^L e -i{**^°)* W here 7:=^. (3.22) 
|Ren 4 Ren 4 
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Figure 4: Graph showing the required negativity of the quantity Re 77,4 as a function of 
h for various values of /. Here we have taken L — 1, K — 2/3 and S = 4/15. Other 
geometries were verified as well. 

We have verified the condition Re 77,4 < numerically for a wide range of parameter 
values. See Figure |U 

Summarizing the analysis above, we have shown: 

Theorem 3.2. Fix a choice of parameters K, L,S,h and I such that the first eigenvalue 
of C in (|3.4p satisfies ImAi and 774 given by (|3.2ip satisfies Re 774 < 0. Then taking 
r = Re Ai + e in (|2.ip . there exists a value Eq > such that for all positive e < eo, the 
system (|2.1[) - (|2.7p undergoes a supercritical H op f bifurcation to aperiodic state (ip e ,(j) e ). 
Applying (|3. 16|) to this solution, we see that 



with a e given by (|3.22|) . 

Remark 3.3. Though we do not present the analysis here, one can show that in fact 
this periodic solution is asymptotically attracting. Since there is no unstable subspace as- 
sociated with L\, nearby points off of the center manifold A4 S are exponentially attracted 
to M £ . Then one argues that nearby non-PT symmetric points on M £ are attracted 
exponentially to the PT symmetric part of Ai £ . This type of argument was carried out in 
\22$ and no doubt the same type of argument will work here though we have not checked 
the details. 

Remark 3.4. One can also carry out the bifurcation analysis in the regime where X\ is 
real, which in particular would correspond to parameter regimes where I is sufficiently 
small. As our primary goal in this article is to address issues related to periodic phenom- 
ena raised in J3jj, J5jj and we did not pursue it here. In this case the center subspace 




(3.23) 
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Figure 5: Stationary magnetic vortices. The curves represent level sets of the phase of 
the first eigenfunction U\. The geometry is L = 1,K = 2/3,6 = 4/15. In (a) we used 
h = 15, while in (b) we used h = 17. In both cases we took / = 10 which is below 
the critical value I c , so that a stationary vortex solution emerges to the full problem 
(I2.1[) - (|2.7I) with leading profile given by u\, cf. Remark I3T41 

is simply spanned by u\ and a stable stationary state emerges for e > of the form 

V> £ ~ Ce 1/2 Ul (3.24) 

for some computable constant C . In the case where h is positive and sufficiently large, 
while Ai is real, one expects this stationary state to have magnetic vortices. Such a result 
might be compared to the single vortex stationary solution found in [91 for a similar model. 

4 Vortex formation 

As we mentioned in the Introduction, vortices form in this problem due to two separate 
effects. One type of vortex, that we term a magnetic vortex, is well-known. Magnetic 
vortices form as a result of the applied magnetic field, and can appear even when 1 = 0. 
We present two examples of such vortices in Figure [5] In both cases we used a rectangle 
with parameters L — 1,K — 2/3,6 = 4/15. In Figure [5^ the applied field is h = 15, 
while in Figure \5jp the applied field is h = 17. In both figures the current is I = 10 
which is below the critical current I c for this geometry and these values of h. Hence, 
the first eigenvalue of the linearized problem is real and so the center manifold here is 
one-dimensional. This is the regime discussed in Remark l3.4l There are other examples, 
not shown, where vortices form even when 1 = 0. 

What makes the present problem unusual is the formation of kinematic vortices, 
that is, vortices that are created even in the absence of magnetic fields. Physically, 
these vortices, just as the magnetic vortices, are points in space-time where the order 
parameter vanishes, the order parameter has a nonzero degree around these points, and 
large phase gradients occur near them. The formation of such kinematic vortices was 
extensively studied by Bendiyorov et al. [6] using numerical simulations of the TDGL 
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equations. The authors report on an unusual effect, where vortices appear periodically 
in pairs along the center line x = of the rectangle, and move along it. They also 
found that, depending on the parameter values in the problem, vortices can either form 
at opposite sides of the boundary and annihilate inside or else nucleate together at an 
internal point on the center line and move away from each other towards the boundary. 
When h = 0, both vortices appear a t the same time and they vanish, either by crossing 
the boundary simultaneously, or by annihilating each other, symmetrically about the line 
y = 0. When h ^ the y-symmetry is broken: the creation of these kinematic vortices 
can take place at different times, and their motion is not symmetric with respect to the 
line y = 0. 

We will use the theory developed in the preceding section to give a simple explanation 
for the formation and motion of kinematic vortices. Our analysis also gives a simpler 
means to compute when and where they form. In fact, we derive an explicit equation of 
motion for the kinematic vortices, fully based on the leading eigenfunction of the operator 
C defined in equation (|3.4I) . After deriving the equation of motion below, we demonstrate 
the different types of vortex creation and motion. One benefit of the new theory is that 
it allows us to easily detect additional types of kinematic vortex patterns, not observed 
in 0. 

Mathematically, the formation and motion of the kinematic vortices are a consequence 
of the PT-symmetry of the problem. To see how they are created and move about, we 
consider the leading order term in the center manifold (cf. (|3.23p ): 

%j) = a e {t) Ul +a £ (t)*u\. 

It is convenient to introduce the notation 

of u 1 (0,y)=g(y)e i ^, 

where £ and \ take the values provided in equation p.22j) . Therefore, along the rectan- 
gle's central line the order parameter is given (at leading order) by 

V(0, y, t) = 2Zg(y) cos {- X t + (3(y)) . (4.1) 

Hence, the order parameter vanishes on the central line x — whenever the equation 

X t = p(y) + tt/2 + mr, n = 0,±l,±2,... (4.2) 

holds. Equation (|4.2[) is the equation of both motion and creation of kinematic vortices. 

Remark 4.1. In the computation performed in this section we replaced the normalization 
condition VS. 6\) with the normalization ui(0,0) ~ 1. This condition has the advantage 
that for all parameters we have (3(0) = 0, and therefore it is graphically easy to compare 
different f3 functions. 

The boundary conditions on u\ imply f3'(±K) = 0. Therefore, three simple shapes for 
f3 might be expected: upward hump, downward hump, or a monotone shape. However, 
our numerical study shows that, while indeed each of these shapes can occur, depending 
on the problem's parameters, other shapes are also present. Also, when h = 0, the 
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(a) (b) 

Figure 6: Different shapes for the function j3(y). The geometry is L — 1,K = 2/3,6 = 
4/15. In (a) we set h — 0, while in (b) we set h = 0.05. 

problem is symmetric with respect to the y axis, and therefore /3 is either an even function, 
as in Figure [6^, or an odd function, if it is monotone. On the other hand, when the 
magnetic field is turned on, and h ^ 0, the symmetry of (3 is broken. 

To demonstrate some of the different possible shapes of f3(y) and their dependence 
on the parameters in the problem, we refer to Figure [5] In all parts of this figure we took 
L = 1, K = 2/3, 5 = 4/15. In Figure[6^ we set h — 0. Here f3 appears as an even function 
of y. Note that it changes its concavity depending on the values of the given current 
I. In Figure EJd we observe symmetry breaking in the shape of /3 under the influence of 
a nonzero applied field h ~ 0.05. In particular we point out that for some I levels the 
function f3(y) becomes a monotone function. We are still pursuing an explanation for 
the dramatic change in f3 as / varies from 100 to 105. 

We now analyze four examples to illustrate the effect of the shape of the function (3 
on vortex creation and motion. 

Case 1: We refer to Figure [71 We use the same geometric parameters as before, namely, 
L = 1,K = 2/3,6 = 4/15. We also use I = 25, /i = 0.05. The function (3 is an 
asymmetric downward hump, as depicted in Figure [7^. The first solution of equation 
(|4.2j) occurs when t is large enough so that \t reaches (3(K) + tt/2 and a vortex emerges 
from the boundary. Then, as t increases the solution moves to lower values of y. When 
t is sufficiently large so that \t = P(—K) + ir/2, a second vortex forms at the lower end 
y = K. As t grows further, both vortices move towards each other. When t reaches the 
level where %i = l3(y m ) + ir/2, where y m is the location of the maximal point of ft, the 
vortices collide and they annihilate. This is an example of annihilation of a vortex/anti- 
vortex pair, where we use 'anti- vortex' to refer to a vortex of negative degree. Then, there 
is a time interval where equation (|4.2[) does not hold for any y, and therefore there are no 
kinematic vortices at those times. This scenario repeats itself when %t = f3(—K) + 3tt/2 
and so on. The creation, motion and annihilation of vortices are shown in Figure [TJd for 
a half of a single period. 

Case 2: Using the same parameters as in the previous example, except increasing the 
current / to take the value / = 110, gives rise to a different shape for (3(y). As depicted in 
Figure[5K, it is now a distorted U shape. Therefore, denoting the location of the minimum 
of /3 by y m , when t reaches the value where xt = fiillm) + t/2, a vortex/ anti- vortex pair 
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(a) (b) 

Figure 7: Creation and motion of kinematic vortices, (a) The function /3(y) for the 
parameters L = 1, K = 2/3,5 = 4/15, h = 0.05,/ = 25. (b) The circles describe the 
location of the vortices in the (y,t) plane. 




-0.8 -0.6 -0.4 -0.2 0.2 0.4 0.6 0.S 02 14 0.215 0.216 0.217 0.218 0.219 

y t 
(a) (b) 

Figure 8: Creation and motion of kinematic vortices, (a) The function j3(y) for the 
parameters L = 1, K = 2/3,5 = 4/15, h = 0.05,/ = 110. (b) The circles describe the 
location of the vortices in the (y, t) plane. 



is created inside the sample. As t increases, equation (|4.2j) is satisfied at two locations, 
until a point of time where yt — + 7r/2. After that, only one vortex remains in the 

rectangle, and this vortex eventually leaves the domain when %t — (3(—K) + tt/2. The 
motion of these kinematic vortices is depicted in Figure [SJd. 

Case 3: For the third example we maintain the same parameters as in the second 
example above, except that we increase h to take the value h = 0.2. For this choice of 
parameters f3(y) is a monotone function as depicted in Figure [9^. Now, the vortex is first 
created when \t = P(K) + tt/2. It then travels to the lower end of the center line until 
\t = /3(—K) + tt/2. This motion is depicted in Figure |9t>. 

Case 4: In the fourth example we present a case where /3(y) has both a local maximum 
and a local minimum. The function /3(y) for the parameters h — 20, / = 25 is shown in 
Figure [TOk. This is the parameter choice of Figure [3] as well. The geometry is the same 
as in the preceding examples. Following the scenarios above, if we look along the center 
line x = we expect to see a first vortex emerging at y = K at a time that we denote 
t\. Then, after this first vortex appears, a pair of vortices (of the same degree) appear 
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Figure 9: Creation and motion of kinematic vortices, (a) The function /3{y) for the 
parameters L = l,K = 2/3,6 = 4/15, h = 0.2,1 = 110. (b) The circles describe the 
location of the vortices along the central line x = in the (y,t) plane. 




(a) (b) 

Figure 10: Creation and motion of kinematic vortices, (a) The function /3(y) for the 
parameters L = 1,K = 2/3, 6 = 4/15, h = 20,1 = 25. (b) The circles describe the 
location of the vortices in the (y, t) plane. 
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at a later time, say t 2 , where = P(y) + tt/2 and y is the location of the interior local 
minimum of /3. They move away from each other, until the one moving upward collides 
at a time £3 with the first vortex moving downward. Finally, the remaining vortex that 
moves downward reaches the boundary y = —K and exits the rectangle. This vortex 
creation and moti on is indeed verified in Figure [TUb. 

This example, however, has several peculiar features. In Cases 1 and 2 vortices formed 
or disappeared in pairs of vortex-antivortex structure on the center line. In Case 4, on 
the other hand, the picture is different. We refer to a sequence of snapshots in Figure [Til 
In Figure [Tib we observe five vortices, with only one on the center line. The two vortices 
far from the center line barely move throughout the period of the evolution. We view 
these two "sluggish" vortices as magnetic vortices. The vortex on the center line is the 
kinematic vortex that formed at y = K at t = t\ as explained above. The two vortices 
located on either side of the center line, and near it, which are of the same degree, move 
towards each other. Eventually they meet at t = t 2 on the center line, giving rise to the 
two kinematic vortices that were discussed above, and are shown in Figure [TDb as well as 
Figure [TTb. These two vortices move as descr ibed above until t = £3. Then, the middle 
vortex on the center line meets the upper vortex on the center line. This is shown also 
in Figure [TTb . Then, as shown in Figure [TTH . this pair of vortices split away from the 
center line. This new pair of vortices moves away from the center line and upward, and 
the entire process repeats itself periodically. 

The picture we just outlined indicates that kinematic vortices can move away from 
the center line. We therefore also term such vortices kinematic, namely those that spend 
part of a period on and part of a period off the center line, since their formation and 
motion follow directly from the PT symmetry and the structure it imposes on the center 
manifold. 

5 Discussion 

All of the examples from the previous section illustrate that, at least near the normal 
state, there is a dichotomy in vortex behavior when both applied currents and applied 
magnetic fields are present. The expansion based on center manifold reduction gives 
a partial explanation for this phenomenon, with the kinematic vortices arising in part 
due to the PT symmetry of the problem. In any event, it is clear that the variety of 
possible vortex behavior in this system is far more extensive than that seen in models 
capturing only magnetic effects. In particular, the motion law (|4.2[) . based on small 
amplitude asymptotics rather than large Ginzburg-Landau parameter asymptotics as is 
more common in the literature, allows for a wide range of effects including boundary 
and interior nucleation, collision of like-signed vortices and periodicity of these events. 
Of course, all of the rigorous analysis we conduct necessarily involves small amplitude 
solutions since it is based on a bifurcation from the normal state. One would imagine 
that an even richer array of vortex behavior is possible for this system if one looks far 
from the normal state though an analytical approach would clearly require different tools. 
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Figure 11: Contour plots of the phase of the order parameter for the parameters of Case 
4 at four distinct times during a period of the motion. In (a) we see five vortices, one 
of them at the upper end of the center line. In (b) two vortices that were earlier off the 
center line move into it and meet there. Then they separate and start moving away from 
each other along the center line in (c), with the upper one moving upward to approach a 
south-moving vortex along the center line. In (d) these two vortices seem to collide and 
then veer away from each other. In these figures, the darkest lines are not significant in 
that they only represent a —it to w jump in the phase. The tips of these lines, however, 
represent the vortices. 
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